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By performing accurate ab-initio density functional theory calculations, we study the role of 4/ 
electrons in stabilizing the magnetic- field- induced ferroelectric state of DyFeOs. We confirm that the 
' ferroelectric polarization is driven by an exchange-strictive mechanism, working between adjacent 

, spin-polarized Fe and Dy layers, as suggested by Y. Tokunaga [Phys. Rev. Lett, 101, 097205 

(2008)]. A careful electronic structure analysis suggests that coupling between Dy and Fe spin 
^ ' sublattices is mediated by Dy-d and 0-2p hybridization. Our results are robust with respect to the 

different computational schemes used for d and / localized states, such as the DFT-|-i7 method, the 
Ifeyd-Scuseria-Ernzerhof (HSE) hybrid functional and the GW approach. Our findings indicate that 
the interaction between the / and d sublattice might be used to tailor ferroelectric and magnetic 
properties of multiferroic compounds. 
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cy^ , I. INTRODUCTION 

Multiferroics (MFs) are compounds where long-range magnetic and dipolar orders coexist There is plenty of 
fascinating physics in these materials, due to the strong entanglement of spin-charge-orbital degrees of freedom^, Q 
and a great potential for technological applications has already been recognized 0-01 ■ The coupling between magnetism 
■ and ferroelectricity can be used for sensing applications, but also for memory devices where data is typically stored 
as magnetic information and read out electrically. Recently, several manganese and iron oxides have been shown to 
possess strong coupling; however, ferroelectricity in these materials is rather weak and only the electrical polarization 
' ^ ' can be switched by a magnetic field (but not viceversa, a limitation for many applications). Within this framework, 
DyFeOs is a very interesting compound [7,, |8| because it shows large ferroelectric polarization combined with a strong 
magnetoelectric coupling. Furthermore, the ferroelectric polarization in DyFeOa is induced by the peculiar magnetic 
structure, obtained upon applying an external magnetic field. 

DyFeOa belongs to the the class of perovkite oxides, such as RMO3 (where R is a rare-earth ion and M is a 
transition metal atom). From the computational point of view, their description poses serious problems. These 
compounds are often (strongly) correlated materials, involving d— and /— electronic charge with significant spatial 
localization. First-principles density functional theory (DFT) calculations in the most commonly applied local density 
approximation (LDA) or generalized gradient approximation (GGA) have to face well known deficiencies: the non- 
locality of the screened exchange interaction is not well taken into account and the electrostatic self-interaction is 
, not entirely compensated [l^. Since semilocal functionals tend to delocalize /-states, the /-electrons are often kept 
^-H "frozen" in the core, and the origin of multiferroicity is generally attributed to the spin-charge-orbital degrees of 
^\ freedom of the M sublattice. Although this "standard" approach helped in clarifying many mechanisms leading to 
' ferroelectricity in MFs[l3|, the influence of / electrons on multiferroicity has not been extensively investigated yet 
^ ' by ab-initio calculations, despite several experiments clearly point out /-electrons to play an important role in MF 
properties of RMO3 compounds [Tsl [l6t . This is especially so for DyFeOs, where the electric polarization appears at a 
very low temperature, corres pon ding to ordering of Dy spins [9l-[T]|. Recall that Dy is expected to have a 3-1- oxidation 
^ . state, i.e. a /^ configuration . 
■ - ' In this study, stimulated by the recent experimental demonstration of a magnetic-field-induced ferroelectric (FE) 
phase of DyFeOa (DF0)|9l-[ll|. we performed ab-initio simulations in order to show the key role playe d by / electrons 
in stabilizing spin-driven ferroelectricity. Our calculations support the experimental study of Ref. [l^, Hll and, at 
the same time, quantify the polarization and shed light on the microscopic mechanism (i.e. based on electronic 
structure analysis) of the origin of ferroelectricity. To the best of our knowledge, this is the first ab-initio study of the 
above-mentioned effect. Our simulations were mainly carried out within a DFT-I-J7 approach for localized electrons. 
In addition, we used the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid functional [itI ITsj. which has been shown 
to improve the description of d- and /-electron systems [l^-[2l| over LDA or GGA. Finally, our calculations were 
benchmarked by single shot GW using HSE wave functions {GqWo@H S E); this treatment is expected to give a very 
accurate description of the electronic ground state|2^]. 
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FIG. 1: (Color online) Ball-and-stick model of DFO in the PE and FE states (FEi and FE2). Upper part: side views; bottom 
part; top view along the b axis. Blue (large), red (medium) and black (small) spheres show Fe, Dy, and O atoms, respectively. 
Oeq, Oap refer to equatorial and apical oxygens. Arrows represent spin moments. 



II. COMPUTATIONAL DETAILS 



Calculations were performed using the projector augmented-wave (PAW) method [ilJl with the Perdew-Burke- 
Ernzerhof (PBE) GGA functional!^. We used DFT+U within Dudarev's approach0|using U^b = U - J=3 eV 
and 4 eV for Fe-d and Dy-/ states, respectively. The energy cutoff was set to 400 eV and a 4x2x4 Monkhorst-Pack 
grid of fc-points was used. We treated the Dy / electrons both as valence and as core states jl^l- The Berry phase 
approach[2^ [2^ was used to calculate the macroscopic polarization P. Non coUinear magnetism was treated in 
accordance with Ref. [s^. Spin-orbit coupling (SOC) was included for the end-point states of the adiabatic path, 
i.e. paraelectric (PE) and ferroelectric (FE) DyFeOs, see below. For hybrid functionals, we used the HSE functional, 
recently implemented in VASP[3l]. GW calculations [32, [33| were performed on top of the HSE ground state[22j. The 
experimental lattice constants for orthorhombic DFO were used (space group Pnma, with a=5.596 A, 6=7.629 A, 
c=5.301 A)[3]. Starting from experimental atomic positions, we performed atomic relaxations until residual Hellman- 
Feynman forces were <0.01 eV/A. 

In figure [TJ we show the two FE states of DFO, FEi and FE2, with opposite polarization. These represent a 
simplified version of the complex experimental magnetic-fleld-induced polar states (cfr figure 1 (c) and (d) of Ref. 
[10'), obtained by neglecting the x and z components of Dy and Fe spins. f4l|. In both, the FEi and the FE2 cases, 
the Fe and Dy sublattices show an A- type magnetic structure, i.e. ferromagnetic (FM) and antiferromagnetic (AFM) 
intralayer and interlayer coupling, respectively, with the spins pointing along the in-plane c axis. In FEi, the stacking 
of the spins along the out-of-plane b axis is J,Fe-|Dy-tFe-tDy; the FE2 state is obtained from the FEi state by rotating 
Dy spins by 180°, so that the spin stacking is 4,Fe-tDy-tFe-4,Dy. In the PE state, the atoms are arranged according 
to the Pnma {D^h) space group and the Dy spins are still intra- layer FM coupled, but rotated with respect to the 
Fe spins by 90°. In FEi, the spins of a Fe layer become parallel to the moments on one of the nearest-neighbor 
Dy layers and antiparallel to the other: Dy layers should then displace cooperatively towards Fe layers with parallel 
spins via exchange striction, giving rise to alternating short-long-short-long interlayer Dy-Fe distances. Accordingly, 
a polarization P along the b axis should be generated. In the FE2 state, the flip of Dy spins would cause a reversal 
of P. Finally, in the PE state, each Fe sheet is sandwiched by layers with ± Dy spins, and no interlayer dimerization 
is expected. This microscopic mechanism, which involves frustrated interactions between rare earth and transition 
metal ions and its optimization by exchange striction, was proposed in Refs. [l^, [111. 



III. SWITCHABLE FERROELECTRIC STATES 

Our calculations confirm the above interpretation. First of all, the total energies of FEi and FE2 are degenerate, [i^ 
Furthermore, from the symmetry point-of-view, the rotation of Dy spins from PE to FEi or to FE2 causes a symmetry 
lowering from the non-polar space group 62 {D2h) to a polar space group 33 (C2t,), giving rise to a polarization along 
the b axis. Indeed, the distortions lead to an alternate short-long-short-long interlayer distance, with dpM = 1.898 
A and dAFM = 1.910 A, while dideai = 1.904 A in the unrelaxed non-spin polarized caseHl]. Clearly, the Dy-/ 
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FIG. 2: (Color online) Adiabatic path connecting the PE and the FEi phases through step- wise rotation of the Dy spin direction 
direction from 90 to 0°. See text for details. Units are (eV/cell) for AE, fj,Ccm~^ for FE polarization, A for interlayer distance 
(d). 



electrons play a key role in stabilizing the FE state. In order to prove this, we froze the / electrons in the core and we 
calculated the electronic ground state using the previously relaxed FEi structure. First, the electronic contribution to 
the polarization, P^ie vanishes; then, by allowing the ions to relax, the ionic contribution, Pionic, becomes negligible 
as well, and the final crystal structure is non-polar. In summary, when treating the / electrons as valence states, the 
PE state becomes unstable, the D2/1 point group symmetry is spontaneously broken to and the system evolves 
towards a stable and polar state. If the / electrons are removed from the valence and frozen in the core, the PE state 
remains stable. This unambiguously confirms that / states are a necessary ingredient for ferroelectricity in DyFeOa. 



IV. DISCUSSION 



A. Adiabatic path 



To shed light on the onset of ferroelectricity, we construct an adiabatic path by progressively rotating the Dy spins 
from 90 to degree, i.e. from the PE to FEi state. The results are summarized in figure [5J panel (I) shows the 
energy difference between the FEi and the PE phases, evaluated at the ideal centrosymmetric (CS) ionic structure (in 
blue) and relaxed configuration (in red); in (II), (III), (IV) we show the electronic, ionic and total FE polarization, 
respectively, evaluated at the ideal (blue) and relaxed (red) ionic structure; in (IV) the Point-Charge-Model (PCM) 
estimate of the polarization is also reported (in black); in (V) the change in the interlayer distances along the path (see 
also Fig.[T]) is shown. Panel (I) clearly shows that ferroelectricity in DFO is magnetically induced. In fact, already for 
the CS structure, the spin rotations give rise to an energy gain (blue bars), further enhanced by ionic relaxations (red 
bars). Note that the relaxation of the electronic degrees of freedoms accounts for most of the total stabilization energy 
(cfr blue and red bars). The energy gain increases from left to right in panel (I), i.e. towards a collinear configuration. 
In parallel, Pr^r also increases from left to right, even in the ideal CS structure (blue), as expected for magnetically 
induced MFs|34|: ionic relaxations further increase Peie [see panel (II)]. From panel (III), we see that Pion is opposite 
to Peie and of the same order of magnitude. However, Pion does not fully compensate Peie [see panel (IV)], giving 
rise to a total polarization Ptot of ^ 0.20 /iCcm~^ for collinear spins. The inclusion of SOC confirms Ptot, i-e. 0.18 
^Ccm~^. Notably, this value is in good agreement with the estimate given in Ref. [l^l- HSE also confirms the value 
of FE polarization. Furthermore, note that Ppcm is not only opposite to Ptot, but also smaller in absolute value than 
Ptot- This confirms that electronic degrees of freedom trigger the FE transition. Finally, in panel (V), the evolution 
of the interlayer distances along the polar b axis clearly shows that the Dy and Fe layers are coupled, giving rise to 
spin-driven interlayer dimerization. 



FIG. 3: (Color online) Left: Atomic displacements (arrows) between the CS and the FEi structure. Spins are not shown for 
clarity. Right: Projection of positive DELF(r) isosurface in the ab plane of FEi structure. See text for details. 

B. Ionic relaxations and origin of the ferroelectricity 

In the PE phase Fe, Dy, Oap (apical oxygens), Oeq (equatorial oxygens) occupy the 4a, 4c, 4c, and 8d WyckofF 
positions (WPs), respectively. When the symmetry is lowered to C2«, Dy and Oap change their site symmetry to 4a 
and the Oeq become inequivalent (8d — > 4a + 4a). This is readily explained by considering the local spin configuration 
around OeqS: when the latter are sandwiched by FM coupled Fe and Dy layers, they have two tFe and two fDy 
atoms as nearest neighbors (these will be labelled 0|^^); when sandwiched by Fe and Dy layers AFM coupled, they 
have two ^¥e and two |Dy atoms as nearest neighbors (labelled as Ol^'^). First of all, all OeqS carry a spin- induced 
moment parallel to the neighboring Fe atom. Furthermore, by imposing the FEi spin configuration on top of the CS 
ionic structure, 0\'^ and O^'^ become inequivalent: the 0|^^ has ±0.194 and the 0\'^ has ±0.207 as induced 
spin moment. To rule out any numerical artifact on this small difference, we impose the PE spin configuration on 
top of the CS ionic structure. In this case, all OeqS carry induced spin moments of exactly the same magnitude. 
In passing, we note that all oxygens remain equivalent when freezing_the Dy-/ electrons in the core. Furthermore, 
we performed an analysis of the symmetry breaking distortions [SSl |36|. The mode decomposition confirms that a 
FE mode is involved, called GMA—. In Fig. [3] we show the pattern of atomic displacements with respect to the CS 
structure. The inequivalency of Oeq is subtle: 0|^^ (Olq^) move in such a way to decrease {increase) the distance to 
its neighbor Dy atom. For 0|^^, doy-o is 2.478 A; for 0|'-^, dcy-o is 2.496 A (the same distance in the PE phase is 
2.487 A), suggesting that a weak bonding interaction is active between the FM layers, leading to a polarization along 
b. A useful tool for studying tiny differences in bonding interaction in solid state systems is the electron localization 
function (ELF)[33, HI]. The ELF values lie by definition between zero and one. Values are close to 1, if in the 
vicinity of one electron no other electron with the same spin may be found, for instance as occurs in bonding pairs. 
Here, we consider the difference in ELF (DELF) between the situation when / electrons are in the valence and when 
they are frozen in the core, for the same ionic configuration, i.e. DELF(r)=ELF/^,^, (r)-ELF/^^^^(r). The physical 
interpretation is as follows: positive values of DELF show up in regions where the electron localization is higher, 
i.e. the bonding between FM layers is strengthened. In Fig. [3] (right part) we show a positive isosurface of DELF 
projected into the ab plane. Clearly, it is mainly localized between FM layers and, more specifically, in the region 
between Dy and O^'^ . This points to a bonding interaction between FM layers mediated by 0|^'''. 

C. Electronic structure fingerprint 

A careful inspection of the orbital-decomposed magnetic moments reveals that the Dy-bd states are polarized only 
if the 4/ states are in the valence: the 4/ states couple locally to the 5d spin moments by intra-atomic Af-bd exchange 
interaction [39]. The 5d states are much more extended than 4/ electrons, suggesting that the glue that finally couples 
FM Fe and Dy layers is the interatomic interaction between the Fe 3d states and the Dy 5d states mediated by 
the oxygen p states. In Fig. U] we show the local Density of States (DOS) at 0|^^ atoms and the Dy-d and Dy-/ 
states. In panel (I), a clear interaction between d and / states of Dy and oxygen states is visible (cfr dotted ellipse). 
Note that the same interaction involves the Fe 3c? states as well (not shown in Fig. 0]). Remarkably, the interaction 
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FIG. 4: (Color online) DOS for oxygen and Dy atoms for / as valence [panel (I)] or as core states [panel (II)]. Panel (III) shows 
the /-DOS as calculated by DFT+U, HSE and GoWoQHSE. Vertical dotted lines refer to the Fermi level. Positive (negative) 
DOS values refer to minority (majority) states. See text for details. 

disappears when freezing the / states in the core (see panel (II)), for which the intra-atomic Dy- f-d and interatomic 
0-p hybridization disappears. Panel (III) shows the Dy-/ DOS calculated using DFT+U, HSE and GoWo@HSE 
calculations. In the (relevant) occupied manifold, the chosen U nicely fits the HSE DOS, which, in turn, is rather 
close to the Go Wo calculations. The corrections beyond a mere DFT+U approach show up in the unoccupied states 
by opening the band gap; however, this does not change our conclusions, as far as the polarization is concerned. 

The fact that the Dy and Fe interaction is mediated by O^'^-sp states suggests possible routes to tailor the ferro- 
electric polarization. For instance, compressive or tensile strain along the b axis should increase or decrease the tilting 
of the octahedra, favoring or disfavoring the interaction via the intermediate Ol^"^ states, i.e. enhancing or reducing 
the ferroelectric polarization. 

V. CONCLUSIONS 

Several results emerge from our study: i) the FE state of DFO is driven by exchange striction, confirming the 
qualitative explanation given in RefsllO. I 111 : ii) two degenerate and switchable polar states exist characterized by 
a sign-reversal of the FE polarization (±P) and linked by an adiabatic path (connecting the two ferroelectric (FE) 
states FEi and FE2 to the same paraelectric reference structure) can be obtained through a relative rotation of the 
direction of Dy spins (with respect to Fe spins); in) the coupling between Dy and Fe spin sublattices is mediated by 
Dy-d and 0-2p states; iv) the estimated FE polarization is in agreement with experiments; (v) by freezing the / states 
in the core instead of relaxing them in the valence, we confirm the crucial role played by / electrons in establishing the 
spin-driven ferroelectricity. More generally, our study suggests that / electrons might play an important role in the 
ferroelectric properties of other RMO3 compounds (where 4/ electrons are often neglected in theoretical calculations); 
(v) last, but not least, our electronic structure analysis suggests possible routes to tailor the ferroelectric polarization, 
owing to the strong Dy-Fe coupling via intermediate equatorial oxygens. Further study is in progress to confirm this 
expectation. 
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